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Abstract 

Magnetic field amplification by a fast dynamo is seen in local box simulations of SN-driven 
ISM turbulence, where the self-consistent emergence of large-scale fields agrees very well with 
its mean-field description. We accordingly derive scaling laws of the turbulent transport coef- 
ficients in dependence of the SN rate, density and rotation. These provide the input for global 
simulations of regular magnetic fields in galaxies within a mean-field MHD framework. Using 
a Kennicutt-Schmidt relation between the star formation (SF) rate and midplane density, we 
can reduce the number of free parameters in our global models. We consequently present 
dynamo models for different rotation curves and radial density distributions. 
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1 Introduction 

>• . . . 

\ Simulations of the ISM in a shearing box domain have shown that turbulence driven by SNe 
■ leads to an amplification of the mean magnetic field. Using the test-field method (Schrinner 
et al. |2005j) , we derived transport coefficients relating the mean electromotive force to the 
mean magnetic field (Gressel 2009). With these we were able to reproduce the time behaviour 
\Q | seen in the simulations. Under conditions found in our own galaxy, and assuming a constant 
circular velocity, a rotation rate <;25Gyr~ 1 was required for the dynamo to work. In order 
to further define the turbulence properties as a function of the star formation rate, rotation 
and gas density, we analysed a comprehensive set of direct simulations. Taking these as an 
^ ' input, we here compute global mean-field maps for a set of different model galaxies. 

X: 

5h ■ 2 The mean-field model 

a • 

Measuring test-field coefficients for a wide set of direct simulations (Gressel et al. I2008L I2011|) 
led to the following scaling relations for the relevant diagonal term in the a tensor, 

= 2kms- 1 ^ (SL) (JL 

for the (downward) turbulent pumping described by the antisymmetric part of the a tensor 

/ \ 0.45 / Q \ 2 / \ 0.3 

a rr h = -a^r = 8kms _1 I — I I — I I — 
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for the turbulent diffusivity 



i 



2 



D. Elstner 



and for the mean vertical outflow velocity 

/ \ 0.4 

u z = 15kms _1 ( — ) . 

\(7o J lkpc 

The relations were derived for SF rates, a, varying from one tenth up to the galactic value 
(Jo = 30 Myr _1 kpc -2 , angular velocities between Hq = 25Gyr~ 1 and IOQq and midplane 
densities from O.lpo up to po — 1 cm" 3 . From the simulations, we moreover found a vertical 
gradient of the turbulent velocity V z u' — 40kms~ 1 kpc~ 1 independent of the star formation 
rate, density and angular velocity. 





Figure 1: Radial profiles of the rotation curve, the star formation rate and the corresponding 
distribution of the dynamo numbers and the turbulent diffusivity for model Bl of Table 1. 

We approximate the vertical profiles for the a tensor by a sin(27rz/h) curve with a scale 
height of h= lkpc. The value of r](z) is chosen to be constant for — 1.5h < z < 1.5h 
and linearly growing with a slope of one third outside this range. For simplicity, we assume a 
constant scale height within our models. We also neglect the anisotropic part of the turbulent 
diffusivity, which seems to be of minor importance for the current models. The rotation curve 
is modelled with a Brandt law 

Q = n a (l + (r/m) 2 )-°- 5 - 

Further we modify the vertical wind above h by adding a radial outward velocity of the same 
size as u z . The wind velocities reach values of 100-200 km/s at z=4kpc, which is an order of 
magnitude higher than in the models of Moss et al. 120101 With these input parameters, we 
solve the induction equation 

0B / -* -* -*\ 

= curl yu x B + aB — 7] cutIBJ 

in a cylindrical domain with r < 15kpc, and of vertical extent — 4kpc < z < 4kpc. Defining 
C a = ahrj^ 1 and Cq = Sloh 2 ?? -1 , we obtain a dynamo number D = C a Cn oc f2 2 ' 5 p~ Q ^<j~ QA . 
The pitch angle, p, can be estimated by tan(p) oc (C a /Cn) ' 5 , scaling as f2 -0 - 25 p~ a05 (7 - 2 . 
These estimates show that stronger SF reduces the dynamo number and increases the pitch 
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angle. It is known that the stationary quadrupole solution of the ct£l dynamo exists only in a 
finite range of the dynamo number. Because the final strength of the regular field also depends 
on the saturation process, this estimate does, however, not provide a prediction for the final 
field strength in dependence of the star formation. Nevertheless, this behaviour still opens the 
possibility for radially extended regular magnetic fields. This is because, in an exponential 
disc, SF decays much faster with radius than the angular velocity, and hence the dynamo 
number may be nearly constant over a large radial range. Applying a Kenicutt Schmidt-type 





Figure 2: Total intensity (colour coded) and magnetic field vectors for model A2 and A4 of 
Table 1. 

law, a oc p , we can specify our galaxy model by a radial density profile, which we leave 
constant up to r p and then exponentially decay with a scale length of 3 kpc, as seen in Fig. Q] 
For the nonlinear back-reaction, we use a classical, however anisotropic, a quenching. While 
the usual quenching of only the diagonal terms in the a tensor would lead to solutions with a 
small pitch angle, independently quenching the pumping term can also saturate the dynamo 
by the increasing field advection from the wind. In this case, the pitch angle of the kinematic 
mode can be preserved (Elstner et al. 2009). 

3 Results 



The models rely on a crude approximation of the vertical profiles for the turbulent transport 
coefficients, which still leave some uncertainty in the absolute numbers given in Table 1. 
Nevertheless, the general trend agrees well with the predictions form the local dynamo number 
analysis. The pitch angle, measured from the magnetic vectors (see Fig. [2]) of the polarisation 
map at r p , increases slightly with the star formation rate as predicted by the ratio C a /Cn 
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above. The growth times of the order of 100 Myr tend to increase with the star formation 
rate, but there are exceptions (cf. Al and A2 of Table 1). No large-scale magnetic field 
amplification was observed for a > lOoo in model B3, and in the weakly differentialy rotating 
model A5. Yet, strong starbursts are usually not long-lasting events and therefore the dynamo 
may still operate on longer time scales. The final field strength is not strongly dependent on 
the SF rate, and only the toroidal field is influenced by the difference in turbulent diffusivity. 
The inverse dependence of the dynamo action on the SF activity is mainly due to an enhanced 



Table 1: Model parameters and results 
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Al 
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0.4 


1540 


1.3 


0.25 


12 


12 


0.18 


A2 
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550 


1.1 


0.25 


14 


14 


0.16 


A3 
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2.6 


192 


0.95 


0.25 


17 


17 


0.25 


A4 


2.5 
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86 


0.65 


0.28 


19 


19 


0.35 


A5 
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1 


23 


0.0001 


0.0001 
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Bl 
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1.6 


344 


1.1 


0.26 


11 


14 


0.18 


B2 
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4.2 


121 


0.9 


0.23 


13 


16 


0.49 


B3 
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11.2 


43 


0.0001 


0.0001 
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B4 


2.5 
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54 


0.6 


0.25 


17 


20 


0.47 



turbulent diffusion. This does not necessarily increase the turbulent velocity but may equally 
change the correlation time. In fact, the preference of magnetic arms between the optical arms 
may be based on this very property of ISM turbulence (cf. Rohde et al. I1999|) . 

The main conclusions drawn from the presented set of simulations are as follows: 

• SF rate and rotation determine the dynamo - calling for adequate galaxy evolution models. 

• Low star formation rates favour the dynamo, explaining coherent inter-arm fields. 

• Strong SF may suppress large-scale dynamo action (no vertical fields in the centre). 

• Explaining the radio-FIR relation will require a different type of amplification mechanism 
- at least for the small-scale field. 
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